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We use the 2 particle irreducible Schwinger-Keldysh effective action to set up consistent equations 
for the velocity and pressure correlations of a turbulent flow. We use these equations to derive the 
Kadanoff-Baym equations describing the relaxation to Kolmogorov turbulence in the absence of 
mean velocities. 
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I. INTRODUCTION 

Since the second half of last century field theory techniques have been applied to an ever widening set of 
problems. Turbulence theory has featured prominently in these efforts. In fact, this field has been so active 
that it would be impossible to give here a comprehensive set of references; we refer the reader instead to the 
monographies [l[ and 0, and the review articles UIH- In spite of all this effort, and without disregarding 
the fact that there are diverging opinions about what exactly has been accomplished [5j, |6| , it can be argued 
that the array of tools from field theory which are relevant to turbulence has not been exhausted yet [7[ . 

As a concrete example, let us consider Gioia and Chakravorty's ansatz for the friction factor in turbulent 
pipe flow [8, 9]. Gioia and Chakravorty argue that a global feature of the flow (the friction it produces on the 
walls containing the flow) is directly determined by the spectrum of turbulent energy at small scales within 
the flow. What is at stake is how "hard" (short wavelength) excitations of the system manage to affect the 
"soft" (long wavelength) scales (see also This is a typical problem in the theory of non abelian plasmas, 

such as the one produced after a relativistic heavy ion collision [ill Il2j |. In the high energy context, this 
problem is usually handled by deriving Kadanoff-Baym kinetic equations for the correlations of the relevant 
fields [HI . While not unknown [3] , this kind of app roach has not been extensively developed in turbulence 
theory, except in the theory of wave turbulence [151 ] . 

The Kadanoff-Baym equations themselves are most efficiently derived as an approximation to a more 
detailed set of Schwinger-Dyson equations for the relevant correlations. The Schwinger-Dyson equations in 
turn may be derived from the variation of the two particle irreducible (2PI) Schwinger - Keldysh (or closed 
time-path, CTP) effective action (EA) @. 

The different elements of this approach have been discussed in more or less detail in the turbulence 
literature: the Schwinger-Keldysh technique, or rather the closely related Martin-Siggia-Rose formalism 

nana 

, [T8l | , the use of generating junctionals (whose Legendre transform is the effective action) and diverse 
resummation schemes [l, III Eli Hi! ■ However, they have been rarely applied together. Moreover, we are not 
interested in the 2PI approach as a way to resum a perturbative expansion for the correlations. Our goal is 
to elucidate what can be learnt about fully developed turbulence from the fact that the main equations can 
be derived from a variational principle. 

The subject of this paper is the formulation of the Schwinger-Dyson equations for the mean velocity and 
pressure, and the two-time, two points velocity-velocity, velocity-pressure and pressure-pressure correlations 
for incompressible flow as derived from the 2PI CTP EA For simplicity we shall consider flow in an 
unbounded domain, and we will be mostly concerned with the regime of very large Reynolds number Re. 
As an application we shall derive a transport equation describing the evolution of a nearly homogeneous 
flow (we define below what we mean by nearly homogeneous) towards a Kolmogorov-like cascade. These are 
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highly formal issues, but they need to be disposed of before more physical questions, such as the effect of 
boundaries and/or shear on the flow may be considered. 

This paper is organized as follows. In next Section we present the problem and discuss the notation to be 
used henceforth. 

Section III is devoted to the subject of random Galilean invariance. This symmetry, which to the best of 
our know ledg e has been introduced by Kraichnan [2l( , is the most important structural feature of the model. 
See [2I. lil Il7l [T9I I22I [23l| . We shall show that the self-energies introduced in Section II, which represent all 
nonlinear effects on the dynamics of correlations, are naturally random Galilean invariant (RGI). This holds 
even when the Navier-Stokes equations themselves are not RGI, due do the random external forces. The 
fact that this important result is almost trivial in this formalism is probably the most compelling argument 
that could be offered in its favour. 

In Section IV we discuss the general features of the correlation dynamics in the case where the mean 
velocity and pressure vanish. As an example, we describe the solution to the Schwinger-Dyson equations 
corresponding to Kolmogorov scaling. 

Finally, in Section V we introduced nearly homogeneous flows, derive the local transport equation which 
describes their dynamics and use it to discuss the approach towards Kolmogorov scaling. This is probably the 
only truly new result in this paper. We present it as a simple application to better judge all the machinery 
introduced so far. We conclude with some brief final remarks. 

We have added two appendices: Appendix A discusses an approximation made in Section VI, and Appendix 
B shows how the Karman-Howarth equation |5j, l24[ and the Kolmogorov 4/5 Law [5j, [25j may be derived 
from the 2PI CTP EA formalism. It seems a fair demand on a new formalism that the main successes of the 
old one should be reproduced without uncommon difficulty. 

II. DEFINITIONS AND NOTATIONS 

In this Section we define the model, introduce the 2PI CTP EA, show that only a subset of the solutions 
to the Schwinger-Dyson equations represent physical flows and, most importantly, settle on the notation to 
be used in what follows. 



A. Navier-Stokes equations for incompressible fluids 

We consider a flow defined on unbounded d-dimensional space. We shall work in eulerian coordinates 
throughout. At each point x and time t the flow is defined by the components U p [x, t] of the velocity and 
the pressure P [x, t] . 

The fluid is subject to an external random force f p [x, t]. It is one of the basic assumptions of turbulence 
theory that at high enough Reynolds number Re the main features of the flow are robust with respect to 
the forcing. We therefore have considerable freedom regarding the choice of the stochastic process f p [x, t]. 
We shall assume it is Gaussian, with zero mean and self correlation 

(f p [x, t] i q [y, t y }) = cN pq [x, t; y, ty] (1) 

The constant c = L d V 2 , where L and V are macroscopic length and velocity scales, is introduced for later 
purposes; observe that c has dimensions. 

There is an important class of forcings, those where 

N pq [x, t; y, ty] = Vf% [x - y] S (t - t y ) (2) 

These make the Navier-Stokes equations random Galilean invariant (RGI) (see next Section). We shall not 
assume the noise self-correlation Q] has the form [2] unless explicitly noted (this will happen in Section V). 
The fluid is incompressible, meaning that 
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V P IP = 



(3) 



The dynamics of the flow is described by the incompressible Navier-Stokes equations 



IP + V 9 [IPIF] + V P P = f p 



(4) 



where Vb is the kinematic viscosity. 

The natural notation we have used so far is too involved for the discussion to come. We shall compress it 
as follows. We first compress the space and time dependence into a single continuous index a 

P [x, t] m> P a (5) 
For the velocity, we replace both the space-time dependence and component index p into a single index j 



IP [x, t] i — ► IP 



The equations of motion [3] and 0] take the form 



(6) 



D, [U, P] = D jk \J k + L ]b P b + g jkl V k U l = f, 



(7) 



L jb W = 

Here j = (x, t,p), k = (y, t y . q), I — (z, t z ,r) and b = (y, t y ). Further 



(8) 



5(t-t y )5(yi-y) 



L 3 b = V p 5 (t - t y ) S (x - y) 



(9) 
(10) 



9jki = 2 i s pq^r + V v «] 5 (* _ h) S ( x - y) S ( f - *«) <5 (x - z) 
All expressions are summed (integrated) over repeated discrete (continuous) indexes. 



(11) 



B. The CTP generating functional 



Following the time-honored procedure we decompose all fields in mean field and fluctuations IP = u J + u J 
and P a = p a +p a . We wish to find self-consistent equations for the mean fields and the two-time correlators 



G 3k = (u j u k ) 



(12) 
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H ja = (u>p a ) (13) 

I ab = (p a p b ) (14) 
The idea is to write both men fields and correlations as derivatives of a generating functional 



sw 

Jjj 



(15) 



etc., where 



e iW/c = ^iSj/c^ (16) 

The brackets denote an ensemble average, and Sj is the source action 

Sj = J+IP + J+P a + h<++W\J k + Kf+U j P a + ^K++P a P b (17) 

By adding the constant c from[I] the new source has the same dimensions as the physical random source 
ij. More explicitly 



e iW/c 



J DWVDP V [f] S (Dj [U, P] - ij) 5 (l^w) e lSj/c (18) 



where V is the probability density functional for the random source f . We shall call the fields in [18] the 
physical fields and denote them by a + superindex. We introduce mirror fields U J_ and P a ~ to exponentiate 
the delta functions. The mirror fields are defined to have the same dimensions as the corresponding physical 
field. By symmetry, we add sources to Sj coupled to mirror fields, products of two mirror fields, and products 
of one physical and one mirror field. 

The result of this activity is a theory of field doublets U J = (U i+ , IP") and P A = (P a +,P a -). Observe 
that the new index J combines the old j of the velocity field and the new index a = ± distinguishing 
between physical and mirror fields, J = (x, t, r, a). Because of Schwinger-Keldysh tradition, we shall call a 
the "branch" index Similarly the new index A combines the space-time and branch labels A — (x, t, a). 

Integrating over the random sources we get 

eiW/c = J DXJDP e i[S+Sj]/c ( 1Q ) 

where 

S = IP- [D ]k \5 k + + g lkl U k +U l +] + W-L ja P a + + W+L ]a P a - + ^IP N \ l L (20) 
N is the noise self-correlation from [TJ 
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C. The CTP effective action 

Our goal is not the generating functional but rather the 2PI effective action V, which is the Legendre 
transform of the generating functional W with respect to all sources. This means that we treat all sources, 
whether they are coupled to a single field or to a binary product, or else to a physical or to a mirror field, 
on the same footing. T takes the form [7[ 

T = S [u J ,p A ] +S 2 - -It In [G] + T Q (21) 
where G denotes the matrix of all 64 propagators 

_ / G JK H JB \ 

^ — I jjAK jAB J 



and 

S 2 = \ {D jk (G k +>i- + &-< k+ ) 

+ 2u j !h u(:' J - + 2g ]kl u k+ (G i+J - + G j -' l+ ) 

+ L ja (H j -' a+ + H j+ > a - + H a+J - + H a - j+ ) + iN jk G k -' j - } (23) 

Formally To is the sum of all 2P1 vacuum graphs with full propagators in the internal lines and a cubic 
vertex gjki [7| . Observe that since only u fields partake in the interaction, Tq is a functional of the G JK 
alone. Also that it is independent of the mean fields. This, and the fact that it is random Galilean invariant 
(RG1), as we shall see in next Section, are its two most important properties. 

Let us summarize what we have said so far. Making a further and most extreme compresion of our 
notation, we may describe our model as a theory of fields with a classical action 

S = i [D M „ + iN^] X»X V + X - 9ilvp X»X v XP (24) 
We enlarge this theory by adding a source action 

Sj = + ±K IU ,X»X V (25) 

and define the generating functional 

e ,ff / £ = / DX^ s+Sj V c (26) 

thereby introducing mean fields 



and propagators 



SW 1 

[x»x u + GH (28) 



5K, IV 2 



G 



The full Legendre transform of the generating functional gives the 2PI effective action 

1 1 C 

T = S[x] + - [D M „ + ilV + g^ pX P] G» u - -Tr In [G] + T Q (29) 

where Yq is independent of the mean fields. 

It is a fundamental property of the Legendre transformation that it may be inverted. This leads to the 
equations of motion 

5r l 

— = [TV + iN M „] x v + - 9lwp [x v x" + G vp ] = -Jp - Kp V x v (30) 

= 2 + iP V + 9^ P x p ~ icG~l - c£^„ } = —Kp V (31) 

where we have introduced the self energies 

dG pu ~ 2 ^ {6Z) 

We see that in this formulation the turbulence model becomes completely analogous to a quantum field 
theory problem [7|, with c playing the role of Planck's constant. In this sense, we may call S the "classical" 
action, as oppossed to the "quantum" action T. Indeed if c i— > the path integral is dominated by its saddle 
points, which means that the evolution of the mean fields is well described by the Navier-Stokes equations 
generated by S. The problem is that fully developed turbulence occurs in the opposite limit c i— ► oo. 

D. Physical solutions 

In the previous sections we have doubled the degrees of freedom of the theory by matching each physical 
field to a mirror field. We have at least two good reasons to do so, namely that in the enlarged theory 
the dynamical equations may be derived from the variation of an action functional, and that the relevant 
symmetry (which is random Galilean invariance) will be easier to work with in the enlarged theory. As we 
shall see in next section, the real part of the S action (the part that is independent of the external forcing) 
is naturally RGI. 

However, only some of the solutions of the enlarged theory really describe physical flows. To see this, let 
us use J for the string of all the sources. Then we have (assuming real sources) 

e -iW[j] 7c = j DVDP exp {_j [5* + Sj ] / c } (33) 
but, writing explicitly the branch indexes, 

- S [\J ]a , P a(} ] * = S [aW a , pP af3 ] (34) 

(no sum over a and 0) so 

- W [J ja , J a /3, Kj a ,kp, K ja ^ aj 3, K aa ^ b p\* = W [—aJj a , —f3J a p, —a(}Kj at kp, -a(3K ja ^ a p, -af3K aa}b/3 } (35) 

Therefore at the physical point J = we must have vt>- = p a ~ = 0. G j+ > k+ , G^- k ~, W+^ a +, H j -' a ~, 
ja+,b+ an( j ja-,6- mugt k e rea i 5 wn ii e Gi +,k ~, H j+ ' a ~, H j ~' a+ and I a +^ must be imaginary. 
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We can be more restrictive about the properties of meaningful solutions. Let us now look more closely 
into the generating functional W [J u ~ , J P ] which is obtained when all sources other than J u ~ and J p ~ 
are set to zero. From the definition [TS] we get 

el w[j u -]/c = J DiDJJD p p [f ] s (p. [ U; P ] _ f . + ju-^ s (lIjTJ* + J Q p -) (36) 

Since the determinants of the operators involved are field independent, we get W [J u -,j p -] = constant. 
This shows that on a physical point we have the stronger constraints 

Q J- = pa- = G 3-,k- = H j-,a- = ja-.b- = Q ^ 

Observe that while u^ - and G^~' k ~ must vanish at a physical point, mixed propagators are meaningful. A 
case in point is the retarded propagator 

Gi+> k - = -ic-j^- = -icG ret (38) 

OJk- 

which will play a leading role in what follows ,26;]. 

III. RANDOM GALILEAN INVARIANCE 

As we have already mentioned, important features of our model may be linked to the way it is affected by 
random Galilean transformations 

U II 03, HH H2, H3- In this section we introduce these transformations 
and investigate the transformation properties of both the "classical" action S and the "quantum" action T. 

A. Classical random Galilean invariance 

We shall now introduce a class of transformations that leaves the real part of the "classical" action I2TJ1 
invariant. Let E p (t) be a time dependent homogeneous vector field, and 

y*(t)= f drW{r) (39) 

Now define the transformation 

U p+ (x,i) i ► W (t) + U p+ (x - y (i) , t) 

/CO 
dr E 9 (r) 9{t-r) V 8 LP+ (x, t) + O (E 2 ) (40) 
- OO 

U*-(x,t) - lp-(x-y(t),t) 

/OO 
dr E g (r) 6(t-T) V 9 U r " (x, t) + O (E 2 ) (41) 
- CO 

P + (x, t) ^ -W (t) [x p - y p (t)] + P+ (x - y (i) , i) 

/CO 
dr W (r) (t - r) VgP + (x, t) + O (E 2 ) (42) 
-CO 
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P-(x,i) » Ep(t)lP-(x-y(t),t)+P-(x-y(t),t) 

/OO 
dr E« (r) (t - t) V„P- (x, t) + (£ 2 ) (43) 
-oo 



We obtain the transformation rules 



V P UP± (x, t) ^ VpU^ (x - y (t) , t) (44) 
V P P+ (x, t) i * -E p (t) + V p P+ (x - y (t) , t) (45) 
V P P- (x, t) i * V p [E,U«- (x - y (t) , t)] + V P P- (x - y (t) , t) (46) 



and 



(x, t) ~ E p (f) + |V+ (x - y (t) , t) - E« (t) V,IP+ (x - y (t) , t) (47) 



This transformation turns the action into 



S [U, P] ^ S [U, P] + l - J dxdt dx'dt' 5 E N pq [x, t; x', t'\ IP" (x, t) U 9 " (x', t') 



where 



(48) 



5 E ~N pq [x, t; x', t'] = N pg [x + y (t) , i; x' + y (i') , t'] - N pq [x, t; x', t'} (49) 
Observe that (5bN P9 = when the noise self-correlation is of the form [2] 

B. "Quantum" random galilean invariance 

Now we discuss random galilean invariance of the full "quantum" action T. 

Let us summarize the results from the previous subsection in terms of the most compressed notation 
introduced in Section II. As we have shown, there is a class of transformations 

X» h-» + C£E a + F£ a X v E a + O (E 2 ) (50) 

where the E a are arbitrary parameters, such that the real part of the classical action is invariant and the 
imaginary part preserves its form 

Sc^S c + ^SN^X^X" + O (E 2 ) (51) 
Observe that this implies the identities 
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|d^" + l - 9lxvp X»Xr} [C£ + F&X"] = 



which may be decomposed into 



(52) 



D^C£ = 



(53) 



(54) 



F^ =0 



and 



We define 



We also assume 



N„„C£ = 



F" = 

-'yet u 



(55) 

(56) 
(57) 

(58) 



making the measure of integration also invariant. It is easy to show that the actual random Galilean 
transformations defined above satisfy this condition. 
The source action transforms into 

Sj - Sj + I J„C£ + [J„F; Q + K^CZ] X»+ l - [K, p F? a + K pv F^ a \ X»X» J F« + O (F 2 ) (59) 



The generating functional itself is invariant, and this leads to the identity 



SS 



5E° 



(60) 



where 



(^) = J " C « + + ^ + ^ + + 2<5N ^ q] ^ + (61) 



We further eliminate the sources in terms of derivatives of the effective action 



K F p 4- -8N 

1 ^p.p 1 va ' 2 t tva 



(62) 
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Finally we obtain the Zinn- Justin identity @, H3| 



If we use the explicit form [29] and the known transformation properties of the classical action, we get 

= -\g» ap G°» [C£ + F£ a x v ] - [D w + iN w + 9fi(7p x° - icG^ - cS w ] F£ a G^ + iS^G^ (64) 
Using the tracelessness of F? a this reduces to 

= -\g^ P G° p [C£ + F^ a x v \ - [D pp + <w P z ff - cE MP ] F&G"" (65) 

Now use [54] 

= ig^pG CTp F^x^ + [g^ p r - c£ pp ] F^ a G^ u (66) 

and [55] 

= cZvpF&G^ (67) 

This equation shows that the quantum part of the action is RGI. When the noise self correlation is of the 
form [2] the full quantum action T is RGI. 

C. Random galilean invariance and self-energies 

We have shown that the quantum part of the effective action is random galilean invariant. This suggests 
it can be written as a functional of the random galilean invariant part of the two point functions only. 

Observe that a random galilean invariant kernel F [x, t x ; y, t y ] must be of the form / (x — y) 5 (t x —t y ). 
Therefore we may regard 

F Ktt [x, t x ; y, ty] = S (x - y - z) 6 {t x -t)S (t y - t) (68) 

as a basis of random galilean invariant kernels. We now seek an inner product where this basis is orthonormal, 
namely 

(F.,t,F*,«) =6(z-z')6(t-t') (69) 

The simplest solution is 

(F G) = J dxdt x dydt y F [x, t x ;y, t y ]* G [x, t } (70) 

where L d is the available volume and we have defined S (t - t'f =Sl5(t- t'). Now we can define a projector 
on the space of random galilean invariant kernels 
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G i ^ G flG7 (71) 

where 

Grgi [G] = J dzdt F ait (F KtU G) (72) 
If Tq is really a functional of Grgi rather than simply G, then it follows that 



[x, t x ; y, t y ] = E rQjS/ 3 [x - y, t x ] 5 (t x - t v ) (73) 



In the particular case of a homogeneous solution, we may drop the t x dependence of Yi ra s p. 

We see that random Galilean invariance puts very heavy constraints on the structure of the self-energies. In 
Section IV we shall use these constraints to investigate the velocity correlations in fully developed turbulence. 

IV. THE STRUCTURE OF THE 2PI EQUATIONS 

From this point on, we shall investigate solutions where all background fields vanish. Therefore, we may 
disregard the background field equations [301 except for the constraints 

V 9 G^+ [x, t; y, t y ] \ x=y t=ty = V,G^+ [x, i; y, t y ] | x=y ^ = (74) 

The dynamics of the model is determined by the Schwinger-Dyson equations 13 1 1 where moreover we put all 
sources equal to zero. 

Our strategy relies on the observation that we have two different ways of computing the inverse correlations 
[G -1 ]^. On one hand, the equations of motion imply 

[G- 1 ] M „ = [D„ w + irV - c£ M „] (75) 

If we split the velocity and pressure fields, this means 



V c J V L A k 

where 



[g-'L = ( - ) Lj n " ) m 



v 1K = D „ + iN JK - cz JK = ^ Dj _-fi + e %_ t _ ) <77 » 

On the other hand, we have both 



and 



r G -n G u P = S p 



(78) 
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G^[G- 1 ], /M = ^ (79) 
From the representation [22] we immediately find 



L AK H KB = Si (80) 



Writing this in full with the help of [TO] we get 



V p H^ [x, t; y, t y ] = -ica a H (x - y) 6 (t - t y ) (81) 
where a is the first Pauli matrix 

-(11) («o 

which here plays the role of a "metric tensor" for the branch indexes. The point of [81] is that it has no 
dynamical information, though the solution may depend on the boundary conditions on the flow. For flows 
in unbound space, as we are considering here, we may write down the solution right away 

H** [x,t;y,y =J**L jH*-*')-»{*-f)]H^ [kM (83) 

H pa <P[k,u] = -c%a a/3 (84) 
k 2 

similarly we get 

V p G^[x,i;y,y =0 (85) 

which means that all possible 36 velocity correlators are transverse. 
The velocity correlations are determined by the equations 

V JK G KL = icAj (86) 

where 

A L = 8 L j + -L,, B H BL ee 5i A 1 . (87) 



c 



This operator admits a Fourier transform analogous to I83[ with kernel 



A" =5 q - 

P P f~2 

For a physical solution, we must have G p ~ ,q ~ = 0. Therefore 



Since the second equation says that G fe+ ' ! is regular, the first implies that T,j +ik+ 
The final two equations now read 

[D j+tk _ -cS J+ , fe _]G fe -^ + =icA l j 



[Dj_ tk+ - cEj-, k+ ] G k+ ' l+ + \iN jk - c£j-, k -] G k ~' l+ = 
To extract the meaning of these equations, let us write 

Vp^qf) [x, t; y, t y ] = V pa , q p + V£4,, Qi/ j + V v q B p . aJj + V x p V v q C a 

where 

^qT^pa,qf3 = ^pT^pa,qf3 = ^ V q A q ,u,{i = Vp-Bp )Ce ,/3 = 

Then the above equations decompose into 



P. ./. • G k = icA^- 



which admits the formal solution 



or else, defining 



T>j-, k+ G k+ > l+ +T> j - tk -G k -' l+ = 



G k+,i+ = l _ G k+,j-v G k-,i+ 
c ' 



G k ^~ = -icG& 



rik-,l+ _ _j r .r<ki 
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Vj-, h - = iM 3k (102) 

the physical velocity correlation 

G k+ - l+ = cG k r i t N jk Gt dv (103) 

which may be regarded as the fluctuation-dissipation theorem for homogeneous turbulence [2(| . The velocity 
fluctuations described by G ++ are equivalent to the stochastic fluctuations of a fluid under a renormalized 
force fj e ™ with self correlation 

(f; e "fT n ) = cMjk (104) 

Analog considerations may be done for the pressure-pressure correlations. In this case the relevant equation 

is 

V JK H KB + L JA I AB = (105) 
linking the pressure-pressure correlations to the longitudinal part of the self-energies. 

A. Fully developed turbulence 

We have already introduced the macroscopic length scale L and velocity scale V, in terms of which we 
have the constant c = (L d V 2 ). We say that turbulence is fully developed when the statistical properties of 
the flow are determined by the scales e = V^L^ 1 and vt, — LV/Ke alone, where Re is Reynolds' number. 
Alternatively, we may introduce the microscopic length and velocity scales 

^V /4 L 



M lfi = (107) 



In terms of these scales we have 



c = L d V 2 = Rc 3rf+2/4 A<V 2 (108) 

The limit Re i— > with A and u* fixed corresponds to laminar flow. In this limit c i— > and the path integral 
is dominated by its stationary points. In the opposite limit choo and there is no simple way of computing 
the path integral. 

In the inertial range, the flow becomes independent of the viscosity iv We shall take as an experimental fact 
that many properties, such as equal time two and three point functions, are "universal" , meaning that they are 
determined by the single dimensionfull quantity e [1[. Moreover, we shall also accept as an experimental fact 
that for two and three point correlations those properties can be derived by simple dimensional analysis. Here 
"simple" means that we disregard anomalous dimensions, which may be important for higher correlations 

Let us assume a traslation invariant, steady flow where all background fields vanish. Then the H — h 
propagators can be written in terms of their Fourier transforms 
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G++(x,t;x',0= / -^|^e l [ k ( x - x ')-(*-'')] G ++[ k ,H (109) 



(2nf (2tt 

The equal time correlations have Fourier transforms 



G++[k] = J |^G++ [k,c] 



(110) 



The dimensions of G+ + are k d V 2 , and since the only velocity scale that can be built out of e and k is 

1/3 

(e/k) we find the Kolmogorov spectrum 

C++ [k] ~ GK£ 2/3 fc^^ (2/3) (111) 

We shall now extend this analysis to all the velocity correlations and the corresponding self-energies. We 
assume the flow is homogeneous and isotropic. This means all correlations can be written in terms of their 
Fourier transforms as in 11091 Moreover we can transversality explicit by writing 

G P^qp = c &P<ig<x,P (n2) 

Let us recall that G pa > qf3 [k,w] has units of V 2 k~ d U3- 1 and therefore Q a ^ [k,oj] ~ cG ra > sl3 [k,w] has units 
of a;" 1 . 

Since the self-energy Ej_,fc+ has the structure [73l its Fourier transform 



E p _,,+ [k, u] = ~ | A pq [v [k] + zfi [k]] + [k] 
It follows that the equation for Q ret = iQ^ , in the limit z^s, i — >■ 0, is (cfr 



(113) 



[i (uj — Q, [k]) — v [k]] Q ret = — 1 (114) 

Recall that Q adv = Qret- 

Since the retarded propagator is real and causal, we must have fi = and v > 0. Moreover v has 
dimensions of frequency, and therefore we must have 

v[k] = Vo (k 2 e) 1/3 (115) 

where vq is a dimensionlcss constant. 

We have already shown that the velocity fluctuations described by G ++ are equivalent to the stochastic 
fluctuations of a fluid under a renormalized force fj era with self correlation Af given in !1021 If the bare noise 
self correlation N has the random galilean invariant form [21 then the same will be true of the renormalized 
noise correlation TV. Projecting over the transverse part, we can write [T03l as 

g ++ [kM = c M[k] 2/ , (lie) 



Direct integration yields 
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r M[k] 

and so the Kolmogorov spectrum implies the well known result 



N [k] = N k~ d (118) 

where 

A^o = — C K e (119) 

c 

V. KADANOFF-BAYM EQUATIONS FOR NEARLY HOMOGENEOUS AND ISOTROPIC 

FLOWS 

We now have all the necessary elements to consider the regression to homogeneity and isotropy in high 
Reynolds number flow. We have in mind a situation where there are no mean flows and the fluctuations are 
already nearly traslation invariant. In such a situation it is convenient to introduced centroid coordinate 



learly trasiatic 

and time @, El Q HI X = (1/2) (x + y) and T = (1/2) (t x + t y ) and write 



G"*<" [x, t x ; y, ty] = J rfk , ^ e i[k(x-y)- w (t*-t 9 )] GP /3, 97 [k) Wj Xj T] (12Q) 



to obtain the partial Fourier transform of [G l \ v from that of G up . Of course, we cannot do this exactly; 
we shall be content to find the adiabatic expansion of [G -1 ]^, meaning an expansion where the different 
terms are classified according to the number of derivatives with respect to the centroid variables. We can also 
derive the inverse propagators from the variation of the 2PI CTP EA. Demanding that both computations 
of the inverse propagators yield the same result becomes a set of equations for the partial Fourier transforms 
of the propagators. Truncating these equations to first adiabatic order yields the so-called Kadanoff-Baym 
equations. 

The pressure-pressure corelations are slaved to the longitudinal part of the inverse velocity-velocity corre- 
lations and we shall not consider them further. We also have the constraint 



dk dtu 
(27r) d (2tt) 



1 d 
p + 2 dXP 



G p+ ' q+ [k,w,X,T] = (121) 



which means that the velocity correlations do not induce a nonzero mean velocity. 

As a matter of fact, the analysis of the first half of Section IV holds, because we did not assume homogeneity 
there. The only new element is that we must derive an expression for the adiabatic expansion of the 
composition of two nearly homogeneous kernels, and we have to see what constraints are necessary to 
enforce the transversality of the velocity correlations, [85] 

Consider two kernels 

F[x,^;y,t,] = J ^.|. e *M-«.-M] f[kii;i X,r] (122) 



and 
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G[ X ,t x ;y,t y ] = f * 7 ^e < W^)-'-(*-- t »)]G[k jWj X,r] (123) 
7 (27r) ( 27r J 

Write 



F * G [x, t x ; y, i y ] = J dzdt z F [x,t x ;z,t z ] G[z,t z ;y,t y ] 



^ ■ ^ e i[k(x-y)- w ( t *-t s )] (j p ^ G) [ k ^ 5 X,r] (124) 



The required formula is [3] 



(F*G)[k,w,X,T] = FG-|{F,G} (125) 



where 



{F,G}=V k F\7 x G- — — -(F~G) (126) 

Let us now consider the transversality constraint. In the strictly translation invariant case, the Fourier 
transform G ra ' sl3 [k, oj] of the propagators can be written as a linear combination of the two tensors A rs and 
k r k s /fc 2 . The transversality condition enforces the vanishing of the latter term. To first adiabatic order, we 
may add a new independent vector as the gradiant with respect of the centroid variable of a scalar function. 
Projecting out the transverse part of such a vector, we may write 



k k z k 



(127) 



where Q" , G^t an d Qzt are quantities of first adiabatic order. The transversality condition 



ik r H — 



1 d 



2 dX r 



G r«,sf) [ k;W5 X,T] = 



(128) 



becomes 



- A rs - 



1 d 



2dXf 



a,0 



•tA' 



sr>a,l3 



(129) 



so 



G^ = 



(130) 



and 



^ 2dX^ (131) 



The second transversality condition 



1 d 

-lK s + -- 



2<9X S 



G ra ' sl3 [k,w,X,T] = 



implies 



g^ = -g^ 

^3r ^2r 



Let us make a corresponding expansion for the inverse propagators 



Then we get 



V pa , q0 = A rs V a ,, + ^V l :j + A [k r A ts 



v 2u -,+g+-- +v.—g+'- - ^-v 



long 



ax* 



If we introduce the retarded propagator 



then 



\v,g+<-] = o 



iQr 



£>-,+ = - [Gret] 



<9X" 



P_ 4. + V 



long 



We are almost done. Folowing [?J we define the density of states 



V = — Im [g ret ] sign (w) 

7T 

For example, Kolmogorov turbulence corresponds to 



V 



1 U 



K 



+ v [kY 



and also 



7 = - Im [-&ret] sign (w) 
which in the case of Kolmogorov turbulence reduces to 

7a- = — 

Next define the distribution function F\ from 

The transport equation is the identity 11031 written as 

Working out the succesive compositions, we get 

N - |^-,+| 2 Q +,+ = y [Di t+ {i>-,+,g + ' + } -v^, + {v*_ + ,g+' + } + 

Next write 



D_+ =R-T 



V*_,=R + T 



So 



Now 



M - (r 2 - r 2 ) g+-+ = -< [r {r, e+<+} - r {r, + {iz, r} e+<+] 



-r. 1 (-*) r / X 



g + - + = t4 ?) ^, FisignM 



(f 2 -r 2 ) 



So finally, if 7^ 



where 



A {JZ, Fi} + F {r, Fi} = [TV - 7r7.F1] sign (w) 
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.4 



(r 2 - r 2 ) 



(154) 



B = 



(i? 2 - r 2 ) 



(155) 



Eq. I153l is quite involved because A, B, M and 7 are themselves functional of F\ through the self-energies. 
However, if we are only interested in the regression to Kolmogorov turbulence, we may approximate them 
by their values in a Kolmogorov-type flow. We show in Appendix A that this approximation is consistent 
with the so-called fourth-cumulant discard approximation Then we get 



R = v (k 2 e 



1/3 



(156) 



r = ibj 



(157) 



A 



U3 



(158) 



B 



So (recall uj = k°) for id > 

ii 2 \i/ 3 
v a oj (Are) 



1/3 



vl (fc 2 e) 2/3 + u< 



(159) 



*/ 2 (fc 2 e) 2/3 W 



5Fi 2 w 
9T + 



(k • V x ) Fi 



(160) 



This equation describes the approach to Kolmogorov turbulence of an initially weakly nonhomogeneous flow. 



VI. FINAL REMARKS 



The goal of this paper has been to bring to the problem of turbulence several tools from nonequilibrium 
quantum held theory Q which, in our view, has not been yet exploited to their full potential. Foremost 
among these is the 2PI CTP EA approach and the possibility of depicting strongly coupled nonequilibrium 
fields by means of Kadanoff-Baym equations [l3| . 

The reason why the method is promising is because the 2PI CTP EA allows to take full advantage of the 
symmetries of the theory, in this case random Galilean invariance. This allows to simplify the problem to 
the point where simple dimensional arguments can be used to completely specify relevant quantities, such 
as the velovity correlations at unequal times. 

On the other hand this paper only sets the stage for the generalization of these methods to the mode 
demanding problems of sheared and bounded flows. We intend to continue our research in this direction. 
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Appendix A: Higher correlations and variations of the self-energies 

The above analysis shows that the self-energies are determined by the velocity and pressure correlations in 
a homogeneous solution. We shall now see how higher correlations may be used to determine the variations 
of the self-energies with respect to the velocity correlations. 

Let us begin with the derivatives of the mean fields 



icG pv 



(161) 



6K vp 2 
and the propagators 



' ' {{X»X v X p ) - x? [x v x p + G v p]} = - {C^ p + x v G»» + x"G^} 



(162) 



xnpv x 

-J17 = sr/ x ' x ^-^ 



ic{(X»X v X p ) - x p \x p x v + GH - x v G pp - x p G vp } = icC pvp 



(163) 



SG" V 
SK, 



pa 



SK, 



[{X»X V ) - x p x u ] 



pa 



IC 



— {(X p X u X p X a ) - [x^x v + G pv ] [x p x° + G pa \ 

-v ^ppa + -a G pp + y.PQP.a^ _ -p ^pa + -a qu P + pj}PQ"<r]} 



IC 



= - {C pvpa + x a C pvp + x p C pua + G pp G vrT + G pG G vp } 



(164) 



We can now take the derivatives of the equations of motion. At zero sources we get 



PV + ilV] G v ° + \ 9tlvp \x v G p ° + x p G v ° + C» pa ] = \ 



(165) 



which yields the explicit representation 



— 2 SpapG 3 P G T l 



(166) 



and 



r -\ r -\ , -^M" 



pr ipv 



5G Tl P 



c: 



Tif<7 



Q pa 



(167) 



\Ppv + *rV + g piJp x p ] Cr T + ^9pu P {Cl pTa + G pT G v ° + G pT G»°} = 



(168) 



icgpupC^ T + 



Cr™ + «T77^ \G™ a G VT + G™ T G V<J ] = 
4 SG^f 



(169) 
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A. Self-energy variations 

We shall now go back to the problem of estimating the self-energy variations in isotropic homogeneous 
turbulence. We have already seen that the assumpton that Tq depends on the propagators only through 
the random galilean invariant component V [G] implies that the self energies can be written as 



£ m , s/ 3 [x, t x ; y, t y ] = J dzdt S m , S/ a [z, t] S (x - y - z) S (t x - t) 5 (t y - t) 



(170) 



where S is itself a functional of V [G] . A second variation yields 



6G"-y te; v,t v ] 



Fra,sl3:u~/,v5 [x — y, £ — T), t x — if] <5 (t x — t y ) 5 (if — t v ) 



where we have already assumed homogeneity in time. 
Our starting point is eq. 11691 rewritten as 



Now use eq. 11671 



r -\ r -\ r <JT p _ ■ r r\p - ^Xip r oTp 



°4 ^Atr^V-r 



to get 



SG X ^ 



g^vp9\^riG pr] + -G fl ^G tp lC 4 v G x ^G^l 



2 

ic 



5E 



SG a 



+ 1 



Act ifjT 



(171) 



(172) 



(173) 



(174) 



Observe that the first and third terms only turn on at k = 0, while the second and third contain the fourth 
cumulant. Therefore, under the fourth-cumulant discard approximation we conclude that the variations of 
the self-cncrgics around a Kolmogorov solution vanish, as we have assumed in the text. 

Appendix B: Three point correlations in homogeneous isotropic flows 

In order to use these equations to find the variations of the self energies with respect to the propagators, we 
need further information about the properties of the higher correlations. In this section we shall investigate 
the three point correlations. 

If all arguments correspond to pressure fields, then translation invariance implies 



a 



where C^' 1 is a scalar with the symmetries appropiate to the coresponding correlation. 
In the case where one index corresponds to a velocity field, 



tx \ y? ty] Z, t z 



dp dujp <iq duj q 
(27r) d (2tt) (2irf W) 
i[p(x-x)—u p (t„-t z )] e i[q(y-z)-w g (t B -tz)]£»a,0>7 



p,w p ;q, uj q \ 



(175) 
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V^C r 3 a ^ p [K,t x ;y,t y ; Z ,t z } = (176) 
This immediately shows that the correlations of one velocity and two pressures must vanish 

C r 3 a ^[x,t x ;y,t y ;z,t z ] = Q (177) 
For two velocities and one pressure let us write 

r r a , s p,~i r . . . . . i _ f dp duj p dq du q 

e i[p(x- Z )-c p (t,-t a )] gibCy--)-^-*,)]^,./^ ^ ^. q w j (1?g) 

Then asking for transversality and that each index corresponds to a vector (as oppossed to a pseudo- vector) 
field we find 

Cp s ^ [ P , q] = A$A$ {4^'^ + S^'Vp"} (179) 
For three velocities (jra,sP,ti we wr j^ e 



(27r) d (2tt) ( 2 ^) d (2tt) 
e i[p(x— )-«,(*.-*.)] e »[q(y— )-«,(t 1( -t,)] C ra,»/3,*7 [ pj Up . q> ^] ( 180 ) 

In this case we have the further constraint 

(p + q) t ^ a,a/3, * 7 [p,q]=0 (181) 

Also the definition of e 

e = V (:c)s C3 r+ ' s+ ' r+ [x, t; y, i; z, f] | x=y=z (182) 

becomes the sum rule 



dp dui p dq dui q ^ r+iS+ir+ 
(2ir) d (2tt) ( 27r ) d (2tt) " 



. ; u,jj u,Li/p u,»4 ULi/q r _|_ s _|_ r _|_ r ,,„,< 



B. The von Karman-Howarth equation 

Let us now consider the equation, which is a particular case of 11651 when x v = 

[T> Ju + iN Jv ] G v ° + \gj vp CT = U'j (184) 

We have two choices for J, J — x, t x ,r, ±, and four for a, a = z, t z , ± or a = z, t z , s, ±. The most 
interesting is J = x, £ x , r, — and <r = z,t z ,t, +. The only subtlety is in computing gjvpC^ pa /2, Observe that 
if J = x, t x , r, — , then necessarily v = x', i x , r', + and p = y, s, +. Therefore we get 
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\ {V (x , )t Cl + -' t+ >° [x', t x ; y, i x ; z, i,] + V ( , )t C^ +;r+ ' CT [x', t s ; y, i x ; z, *,]} | x , =y=x (185) 

Let us Fourier transform these expressions with respect to x — z and t x —t z . If J = x, t x , r, — and a = z,t z ,t,j 
we get 

lks I ^W) c;+ ' s+ " 1[k ~ q,UJ ~ LUq;ci,LUq] (186) 

and finally 

DG r+ > t+ + iN rs G s ~> t+ + iK f ^*^±C r 3 + - s+ ' t+ [k - q, w - q, «,] - (187) 

Let us make the index dependence explicit. We assume the bare noise self correlation is transverse and 
delta-correlated in time 

N rt [k,w] = A rt N b [k] (188) 

From the symmetry of the propagators, 

duj uj G r+ ' t+ = (189) 
In the inertial range, both the bare forcing and the viscosity are negligible. We must conclude that 



/ dq duj duj' r+s+t+ , 



For all k in the inertial range. On the other hand, from the sum rule 11831 



dk f dq duj dw' + + + 

J k s / - — ^T^-TTTT^a [k-q,w;q,w] =e (191) 



2nf 8 J (2n) d (2tt) (2tt) 3 



and so we are led to 



/" dq g?w do; r+ s+ r+ . , 
?^(^)(^) C3 ' [k-q,-;q,-] = (2 7 r) e«5(k) 



(27T) 

From isotropy, symmetry in (r, s) and transversality in z, we must have 



(192) 



Therefore 



A=-(27r) d e <5(k) 



(194) 



C. The Kolmogorov 4/5 Law 



We shall use this equation to verify the Kolmogorov 4/5 Law 



{x.[u(x)-u(0)]} 3 ) = -- £ r 4 



Observe that the left hand side can be written as 



— 6x r x s x i C , 3 +,;s+ '* + [x, t; x, t; 0, t] — — 6x r x,sX t 

= -6x r x s x t 



(2,y 

dk 



( 2 ^) d (2tt) (2tt) 

ikx ( — 



fc 2 



{A rt k s + A st k r } 



= 12 



(2^r 



(a) 



, kx (kx) 



k 2 k 4 



Now let us call 



Then 



and 



1 y (2^)^ ^ 



" ' (2nf k 4 



/ 



(27T) 



fc 4 



r— X -2>(r— V 2r— If 
dr / \ dr ) dr \ 2 



Now 



V 2 F 2 = -Ft 

and 

V'F, = -1 



so 
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F, = ^ (204) 



Which is the desired result. 
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